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SIGNAL PRE-PROCESSING FOR ESTIMATING ATTRIBUTES OF A 

TRANSMISSION T.TNK 

The present invention relates to a method for estimating 
attributes of a transmission line from reflectometry measurements, as 
5 described in the preamble of claim 1, and to an apparatus wherein the 
method is apphed, as described in the preamble of claim 3. 

Such a method is already known in the art, e.g. from the article 
entitled 'Estimation of the Transfer Function of a Subscriber Loop by means 
of a 1-Port Scattering Parameter Measurement at the Central Office’, 
10 pubhshed in the Journal of Selected Areas in Communications (J-SAC) of 
June 2002. Therein, a method to predict the transfer function of a local loop 
from reflectometry measurements is described. 

A local loop is a two-conductor transmission line that connects 
customer premises to a central office of a telecommunication network. 
15 Originally, the local loop has given the customer access to the Plain Old 
Telephony Service (POTS). This copper-wired access network has been de- 
facto the medium to dehver at the best cost new telecommunication services 
to residential users and small business offices. Digital Subscriber Line 
(DSL) services exploits the frequency band above 4 kHz up to several MHz, 
20 which is not used by the POTS, to provide the customer with high speed 
access (up to a few Mbit/s). 

Nevertheless, the local loop was engineered for voice-band 
transmission without any prospect for data telecommunication services 
operating at higher frequencies. The local loop may have several 
25 impairments for DSL services that do not affect the POTS, such as bridged 
taps (open ended line derivation), mixed wire gauge, sphces, etc. DSL 
services are mostly affected by reflections arising along the line, which stem 
from any non-uniformity in the impedance of the transmission line. 
Reflections interfere with the direct DSL signal and, at those frequencies for 
30 which the interference is destructive, reduce the magnitude of the received 
signal or alternatively decrease the Signal Noise Ratio (SNR). It is therefore 
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necessary to evaluate a priori the line quality before implementing DSL 
services. 

In a preliminary step of the known method, a physical model of 
the local loop is defined from the electromagnetic theory specifying how the 
5 voltage and current waves propagate throughout the loop. 

The local loop is modeled as comprising a plurality of line 
segments spliced to each other in a tree-like structure. 

A line segment is modeled by the brace where Zo is the 

characteristic impedance of the fine segment, y is the propagation constant 
10 of the line segment, and 1 is the length of the line segment. Zo and y can be 
expressed as a function of the frequency f and of a fini te number of 
coefficients, the coefficients being defined fi'om the cross-sectional geometric 
dimensions of the line and the material constants of the line. 

specifies the way the voltage and current waves propagate along the line. 
15 The boundary conditions, which the line voltage and the line current must 
comply with at a splice or at a termination, yield the reflection and 
transmission factor of the voltage and current waves at that splice or at that 
termination. Extending the reasoning to the entire loop yields a parametric 
model of the local loop specifying how the voltage and current waves 
20 propagate throughout the loop, such as the scattering matrix or the ABCD 
matrix. 

In a further step of the known method, an initial electrical signal 
is generated and coupled to a port of the local loop, while the electrical 
waveforms at that coupling point are observed. The observed waveform is 
25 comprised of an amount of the initial signed and an amount of each signed, if 
any, that is reflected back from any non-uniformity in the impedance of the 
tremsmission line. This measurement procedure is known in the art as 
reflectometry, more specificedly as time domeiin reflectometry or firequency 
domain reflectometry depending on whether the observed waveform is 
30 measured in the time domain or in the frequency dom ain respectively. The 
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measiired signal is called a reflectogram, if no explicit reference to the 
measurement domain is meant, else a time domain reflectogram or a 
frequency domain reflectogram. 

At that stage of the qualification process, the local loop is 
5 considered as a one-port device, the port at which the reflectometry 
measurements are performed being the reference port. 

In a further step of the known method, the coefficient values of 
the parametric model are estimated. The estimation process includes the 
steps of: 

10 - assigning an initial value to some of the coefficients by making use of a 

quasi-realistic model of the local loop, 

- computing from the parametric model the 1-port scattering parameter iS",, 

of the local loop, namely the 1-port scattering parameter 5,, 

being defined as the ratio in the frequency domain of the backward 

15 travelling wave or reflected wave over the forward travelling wave or 
incident wave, 

- determining from the reflectometry measiurements the 1-port scattering 

parameter .S',, of the local loop, namely , 

- quantifying a difference between and by introducing a cost 

20 function, 

- reducing the difference by adapting the coefficient values of the 
parametric model. 

Both the measured and the estimated 1-port scattering 
parameters are determined in a given reference impedance base. 

25 The quasi-reahstic model provides initial values to the coefficients 

of the parametric model for the computation of the 1-port scattering 
parameter . The initial values of the coefficients should preferably 

have the order of magnitude predicted by their formixla so that the 
optimization of the cost function determines a global minim um and does not 
get trapped in a local minimum. The quasi-reahstic model is a white box 
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model with as much a-priori knowledge as possible, as opposed to a black 
box model. The quasi-re abstic model of the local loop is based upon an 
assumed loop topology, and upon the type of transmission line used in the 
field as presumably known to the telecommunication operator. 

5 Once the coefficients have been estimated, it is possible to 

compute the transfer function between any port of the local loop, the local 
loop now being considered as a multi-port device. 

The known method yields good results up to a certain line length, 
e.g. 2000 meters for 0.4 mm polyethylene cable. For longer loops, the first 
10 direct reflections originating firom the coupling point between the 
measiirement device and the line dominates the measmements and make 
the convergence of the parametric model impossible. 

It is an object of the present invention to ehminate the near-end 
reflections and to make the qualification of longer loops possible. 

15 According to the invention, this object is achieved by the method 

defined in claim 1 and by the apparatus defined in claim 3. 

The major part of the near-end reflections are eliminated by 
carrying out the steps of: 

- estimating the characteristic impedance of the local loop, more 
20 specific£dly estimating the characteristic impedance of the fine segment 

whereto the measurement device is coupled, 

- converting the 1-port scattering parameter of the local loop firom the 
reference impedance base to the estimated characteristic impedance base. 

However, there is still a residual influence of the near-end 
25 reflections. This residual influence is mostly due to the fact that the 
estimated characteristic impedance does not correspond perfectly to the real 
characteristic impedance of the line. For example, the estimation of the 
characteristic impedance does not take into account inhomogeneities of the 
line and the resulting backscatter. 
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It is then highly necessary to determine a time zone wherein the 
residual near-end reflections are enclosed so as their contribution is further 
neutrahzed in the process of qualifying the local loop. 

The processing of the 1-port scattering -parameter in accordance 
5 with the present invention can be performed indistinctly in the time domain 
or in the frequency domain, or in any theoretical referent as known to the 
skilled person. 

The present invention is apphcable to any t5q)e of transmission 
line, irrespective of the natxure of the transmission medium and irrespective 
10 of the telecommunication technology being used thereon. 

The present invention is not restricted to the mentioned prior art 
but is apphcable to any method or apparatus making use of reflectometry 
measvirements to estimate attributes of a transmission hne, no matter what 
those attributes are and no matter the domain in which the estimation is 
15 performed. 

Further characterizing embodiments of the present invention are 
mentioned in the appended claims. 

It is to be noticed that the term ‘comprising’, also used in the 
claims, should not be interpreted as being restricted to the means hsted 
20 thereafter. Thus, the scope of the expression ‘a device comprising means A and 
B’ should not be limited to devices consisting only of components A and B. It 
means that with respect to the present invention, the relevant components of 
the device are A and B. 

Similarly, it is to be noticed that the term ‘coupled’, also used in the 
25 claims, shoidd not be interpreted as being restricted to direct connections only. 
Thus, the scope of the expression ‘a device A coupled to a device B’ should not 
be limited to devices or systems wherein an output of device A is directly 
connected to an input of device B. It means that there exists a path between 
an output of A and an input of B which may be a path including other devices 
30 or means. 
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The above and other objects and features of the invention will 
become more apparent and the invention itself will be best understood by 
referring to the following description of an embodiment taken in conjunction 
with the accompanying drawings wherein: 

5 - fig. 1 represents a measurement device according to the invention, 

- fig. 2 represents an instance of a local loop chosen as a case survey, 

- fig. 3 represents a firont end analog unit of the measurement device, 

- fig. 4 represents a measured input impedance and an estimated 
characteristic impedance for a 5 km and a 1 km long 0.4 mm polyethylene 

10 cable, 

- fig. 5 represents an impulse response of a 3 km long 0.4 mm polyethylene 
cable in 100 Q base and in an estimated characteristic impedance base, 

- fig. 6 represents the removal of the time alias fi*om an impulse response 
by fitting an exponential decaying function on the tail of the impulse 

15 response, 

- fig. 7 represents an equal level impulse response and the determination 
of the neutral zone by means of an extrema and a basehne threshold. 

The test device TEST as depicted in fig. 1 comprises: 

- input means IN, 

20 - output means OUT, 

- a time domain reflectometry unit TDR, 

- a processing unit PROG. 

The input means IN and the output means OUT are coupled to 
the processing unit PROG. The processing unit PROG is coupled to the time 
25 domain reflectometry unit TDR via a data communication bus GOM, such as 
a serial bus. The time domain reflectometry unit TDR is coupled to the local 
loop LP. 

The test device TEST is adapted to estimate the transfer function 
of a local loop from reflectometry measurements. The test device TEST is 
30 integrated as a handheld testing device that can be used by a technician in 
the field, e.g. at the central office next to a distribution frame. 
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The input means IN comprise a keyboard and a computer mouse 
associated with a graphical user interface and are adapted to: 

- trigger the time domain reflectometry measurements, 

- enter or select from a drop-down menu any useful pieces of information 
5 such as characteristics of the local loop (used in the parametric model), 

measurement settings, etc. 

The output means OUT comprise a display screen, e.g. to provide 
a technician with the outcome of the measurement procedure. 

The time domain reflectometry unit TDR comprises the following 
10 functional units: 

- a generation unit GEN, to generate an electrical signal of a given shape 
and within a given spectrum range, 

- a measurement unit MEAS, to measure in the time domain an electrical 
signal, such as a voltage, including an analog to digital converter, 

15 - a coupling unit COUP, to couple the excitation signal from the generation 

unit GEN to a conductor pair and to provide the measurement unit 
MEIAS with the appropriate measurement access points, 

- a low-pass anti-ahas filter FILT for the ansdog to digital conversion, 

- a low-noise amplifier LNA, to amplify the signal from the local loop LP, 

20 - a synchronization unit SYN, to trigger synchronously the generation unit 

GEN and the measurement unit MEIAS. 

The processing unit PROC comprises a Pentium IV processor 
supplied by Intel Corp, as well as all the necessary hardware as known to 
the skilled person, such as random access memory, a hard disk, an 
25 interruption controller, a DMA controller, a graphic card, I/O peripherals, 
etc. 

The processing xmit PROC runs software which carries out any of 
the method in accordance with the present invention, more specificedly: 

- gather and average the measurement results from the time domain 
30 reflectometry unit TDR, 
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- process the measurement results in accordance with the present 
invention, 

- estimate therefrom the transfer function of the local loop, as already 
known to the art. 

5 The electrical signal used to excite the loced loop is a multisine 

wave ranging from 5 kHz up to 1 MHz with a tone spacing of 5 kHz. The 
signal type as well as the burst dviration can be set via the input means IN. 

The measured electrical signals are the line voltage in conjunction 
with the reflected voltage wave. The latter allows the discrimination of the 
10 incident and reflected waves within the line voltage. 

Upon request from the operator via the input means IN, the 
processing unit PROC triggers the time dom ain reflectometry unit TDR via 
the data communication bus COM for it to carry out a time domain 
reflectometry measurement. Thereupon, the synchronization unit SYN 
15 triggers both the generation unit GEN, to generate an excitation signal, and 
the measurement unit MEAS, to start measuring the resulting signals at 
the coupling point. A measurement is a set of digital samples, each sample 
being the encoded magnitude of the measured signal at a particular 
instance of time. The sampling frequency is set to 10 MHz, that is to say the 
20 cut off frequency of the anti-alias filter is about 5 MHz. The capturing 
window is 200 ps long so as the frequency resolution matches the tone 
spacing of the excitation signal. The measurement results are delivered via 
the data communication bus COM back to the processing unit PROC for 
further processing. The measiirement procedure is repeated a few times 
25 over a small period so as to compute the mean and variance value of each 
sample. 

In the following, the operation of the preferred embodiment will 
be described with reference to an instance LP of a local loop. However, the 
present invention is not limited thereto. 
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The topology of the local loop LP, as known to the 
telecommunication operator and as entered via the input means IN, is a 
single transmission line. 

Fig. 2 depicts the electrical circuit of a single transmission line 
from which the parametric model is estabhshed. At one end of the line 
stands a voltage source Vg serially connected to a generator impedance Zg. 
At the other end stands a load with an impedance Zl. The line is 
characterized by its transfer function e~^ and characteristic impedance Zo. 

The hne is considered as a 1-port device, with V and I being the 
input voltage and input current at that port. 

The forward and backwsurd traveUing waves a and b at that port 
are defined with respect to a reference impedance Zref: 

a) 

V-ZA 



b = - 






( 2 ) 



The transmission line t5rpe, as known to the telecommunication 
operator and as entered via the input means IN, is an Unshielded Twisted 
Pair (UTP) consisting of 2 parallel circular cylindrical conductors immersed 
in a homogeneous dielectric, which is assumed to be linear and isotropic. 
(Zo,e"'^) of a UTP is given by: 






(3) 



25 



(4) 

yp^ 

where: 

s — jco — j'lT^ represents the Laplace frequency, 

'P = — denotes an auxiliary function used in 

-t- JfyJ^ 

the calculations, 

Jj = denotes the Bessel function. 
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yp = a^s/l is the parallel per-unit length admittance of the transmission 
hne. 

The coefficients ai to as may be given initial values in accordance 

with: 

5 a 



_ 1 [m g ^ 
' ofT V cr arcosh(D/2a) 




=a^[Jia 



_ 

^ ln(l/V^) 



^5 = ; ; 1 

arcosh(Z)/2a) 

10 where: 

a is the radius of the circular conductors (same radius for both 
conductors), 

D is the distance between the 2 conductor axis, 

cris the conductivity of the conductor material (e.g., copper), 

15 - e and p are respectively the electrical permittivity and magnetical 

permeability of the dielectric material (e.g., polyethylene). 

Prior to the reflectometry measurements of the local loop LP, the 
cahbration procedure takes place. 

A more accurate model of the time domeiin reflectometry unit 
20 TDR is proposed in fi^g. 3. In this model, the following elements are taken 
into consideration: 

- the difference 6g of the actual generator impedance with respect to the 
assumed generator impedance Zg, 

- the coupling unit COUP, modeled as a 4-port coupling device and 
characterized by its own scattering matrix. 



25 
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- the input impedance Za and Zb of the devices measuring the line voltage 
Vline and the reflected voltage wave Vrefl. 

Zin is the input impedance of the device under test. 

The forward and backward travelling, waves a and b can be 
determined from the measvired voltage Vline and Vrefl using the formula: 



'b 








V 




a 






^ 4 . 


Aw 





where A,i, X. 2 , Xs and A.4 coefficients are defined from the scattering matrix of 
the couphng unit COUP. 

The 1-port scattering parameter in the reference impedance base 
can be derived from the formula (5): 



^ measured 





\ ^line ^ ^refl 
line 



/ ^ line 

/ ^ line 




( 6 ) 



The cahbration procedme consist in deter min ing the values of the 
coefficients Ci, C 2 and C 3 . This is achieved by measuring the line voltage 
Vline and the reflected voltage wave Vrefl with 3 distinct values of the input 
impedance Zin for which the resulting 1-port scattering parameter is 
preliminary known: 

- first measurement with an open circuit, Zin = 00 , 5',,| = 1 , 

- second measurement with a short circuit, Zin = 0, 5,,| =-l. 



- third measurement with the reference impedance, Zin = Zref, .Sj,| = 0 . 

20 The processing unit PROC gathers the 3 measurement results from the time 
domain reflectometry unit TDR and obtains therefrom 3 linear equations 
from which the values of Ci, C 2 and C 3 can be derived. 

The processing unit PROC determines the 1 -port scattering 
parameter of the local loop LP in the reference impedance base by using the 
25 formula ( 6 ), the value of the coefficients Ci, C 2 and C 3 as determined during 
the cahbration procedure and the measurement results of the line voltage 
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Vline and reflected voltage wave Vrefl, the local loop LP being the device 
under test. 

The far-end reflections should be made visible as clear as possible. 
They may not be distorted by the mismatch between the reference 
impedance Zref and the characteristic impedance Zo, otherwise very small 
but still very important reflections, which have traveled a few thousand 
meters edong the line, will be drowned in the immediate reflected signal. A 
transformation of the base impedance of Sn is then highly necessary. 

The ideal base impedance is the characteristic impedance Zo of 
the line segment whereto the test device TEST is coupled. Estimating this 
characteristic impedance with respect to the measvired input impedance 
yields good results. Indeed, the longer the line is, the closer the input 
impedance is to the characteristic impedance. 

The input impedance can be derived from the measured 1-port 
scattering peirameter as follows: 

y measured y ^ 

^ j ^measured ^ ^ 



The estimation of the characteristic impedance of the local loop 
LP is based upon the parametric model as expressed in the formula (4). 

The preferred estimator is the Maximum Likelihood Estimator 
(MLE) but the present invention is not limited thereto. MLE is chosen 
because it is consistent, asymptotically not biased, asymptotically normally 
distributed, and asymptotically efficient with realistic noise condition. 

The MLE cost function is given by: 









r measured 



(A)| 



k=\ 






( 8 ) 



where: 

- N is the number of frequency samples, 

- fk is the frequency, 

- 0 is the parameter vector of the parametric model of the characteristic 
impedance: 0 = (ai,a2,a3,a4). 
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- (Jk,d) is the frequency sample of the parametric model of the 
characteristic impedance, 

- is the frequency sample of the measured input impedance, 

■ o’zin CA) is the estimate of the veiriance of the k‘l> frequency sample of the 
measured input impedance. 

The minimum of the cost function as a function of the parameter 
vector 0 is found by means of an iterative algorithm. The preferred 
minimizer is the Levenberg-Marquardt (LM) method but the present 
invention is not hmited thereto. 

The minimizer determines a vector 9 for which the minimum of 
the cost function V is reached. The estimated characteristic impedance is 
then given by Zo . 

As can be seen in fig. 4, the estimated characteristic impedance 
does not take into consideration the higher order variations (peaks and 
valleys around the asymptote), which represents the contribution of the load 
impedance and other fine segments, if any. Those contributions are 
irrelevant as far as the characteristic impedance is concerned and are to be 
over looked by the estimation process. 

The measured 1-port scattering parameter is then converted from 
the reference impedance base to the estimated characteristic impedance 
base as follows; 



^measured 



measured 

^in ~^0 

^meastired ^ 



(9) 



As can be seen in fig. 5, the far-end reflections are noticeable after 
the conversion, yet there are still some non neghgible near-end reflections. 

The determination of the time zone wherein the residual near-end 
reflections are enclosed — named hereafter the neutral zone — , shall be 
carried out in the time domain using an inverse Discrete Fourier Transform 



(DFT) of 5,7 



measured \ 



The resvdting signal is named the impxdse 



response of the local loop. 
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However, inverse DFT cannot be applied directly on 
^measured\ (/) becauso not all frequency samples contain meaningful data 

There are regulations by the telecommimication operators that excitation 
signals below and in the voice band are not allowed on subscriber lines. 
5 Also, excitation of such lines above 1 MHz is not advised due to the extreme 
line attenuation. Last, the measiu*ed signal gets no DC component owing to 
the AC couphng. Subsequently, a number of spectral hues will not contedn 
any imposed power. 

Non-energetic frequency samples of (/) are reset to zero, 

\Zq 

10 while the DC component (q) is set to 1, as predicted from the 

*Zq 

theoretical model. The inverse DFT of such a spectrum will give rise to 
ringing effects in the time domain. The ringing will be smoothed away with 
a hner phase scale-spaced filter or any other appropriate filter as known to 
the skilled person. 

15 The inverse DFT of the filtered (y) contsdns time ahas 

IZq 

because the capturing window T is not long enough to reduce the time ahas 
below the quantisation level. The time ahas wiU be reduce to an acceptable 
level by fitting an exponential decaying function A'e~^‘ on the tail of 

IZq 

(see fig. 6). 

20 A least square estimator is used to find the parameters A’ and B’. 

The parameters are estimated on the last portion of the tail. 

The time ahas g(t) is given by an infini te sum of exponential 
decaying functions: 



g(t) = Ae-"' + Ae-‘‘<‘*^ + + ... = Ae“Y,e 



As 






The relation with the estimated parameters is given by: 



B = B' 
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A = A’e^^ 

The de-aliased impulse response is obtained with: 

( 10 ) 

The skilled person is entitled to consider any other alternative but 
5 the exponential deca 3 ring series and the least square estimator as a vahd 
alternative. 

The determination of the neutral zone is achieved by estimating 
the start of the far-end reflections. 

First of all, the de-aliased impulse response is amphtude- 
10 equahzed to compensate for signal attenuation. 

The equal level impvdse response eqwtievd obtained by 

mxdtip lying the de-ahased impulse response j^auased ^ time 

varying amplification since the more delay, the more attenuation. The 
amplification at time instance t depends on the maximum value of the 
15 impulse response of a generic line of length L = x / , with Vprop being the 

propagation speed of the line. 

The time varying amplification is topped firom a certain time 
onwards to avoid excessive amplification of the tail. 

The skilled person is entitled to consider any other amplification 
20 method but the foregoing time varying amplification, including the non 
amplification of the de-aliased impulse response, as a vahd alternative. 

Next, two threshold levels allow the detection of the far-end 
reflections in equniievei- first threshold allows finding the first 

extremum, which should correspond to the first direct far-end reflection. 
25 Starting fi'om that vahd extremum, the curve of equal levei foUowed in 

the negative time direction. When a crossing occurs with the second 
baseline threshold, the start of that first direct far-end reflection is assumed 
to be found. The neutral zone goes firom the time origin up to that time (see 
fig. 7). 
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The skilled person is entitled to consider any other alternative but 
the foregoing algorithm as a vahd alternative. An algorithm based on the 
computation of a squeure mean value over a sliding window moving onwards 
could be such an alternative. 

5 Ultimately, the so-identified time zone is neutralized in the 

process of estimating the transfer function of the local loop LP. 

Referring to the above mentioned prior art, one could for instance 
set to zero the time samples within said time zone, so as the near-end 
reflections does no longer preclude the parametric model fi*om converging. 

10 Care must be taken to process the estimated 1-port scattering 

parameter in the same way than the measured one was, so as the 
convergence is not biased. 

In an alternative embodiment of the present invention, the 
processing unit PROC would accommodate a feature extraction algorithm 
15 operating in the time domain. The algorithm could use the so-identified time 
zone to distinguish between reflections caused by impedance irregularities 
and reflections coming from line ends, gauge changes and junctions, the 
latter ones only being significant. 

In an alternative embodiment of the present invention, the 
20 processing unit PROC would use a measured curve from a database of Zo 
cxuwes as the base impedance. In order to select the proper curve, the 
technician is requested to enter via the input means IN information about 
the cable type. Yet, the behavior of the characteristic impedance depends 
also on external influences such as temperature and aging. And last but not 
25 least, the properties of a twisted pair vary from pair to pedr in the same 
cable due to the different twisted rates used. Deviations of 50 Q at 10 kHz 
and 5 Q above 100 kHz were measured. 

In an alternative embodiment of the present invention, the 
capturing window would be increased — up to the extent of the available 
30 memory — so as to include the whole impulse response. The foregoing step of 
removing the time ahas would be no longer necesseury. Yet, the resulting 
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frequency resolution may contains samples within some prohibited 
frequency band, maJdng the qualification of some transmission lines 
impossible. 

In an alternative embodiment of the present invention, part or 
5 whole of the processing would be performed on a remote host. For instance, 
the time domain reflectometry unit would be mounted on a card and 
plugged into a ceird slot of an access miiltiplexer, and would send 
measurement results over a data communication network to a remote 
processing unit, e.g. to a Personal Computer (PC) over a Local Area 
10 Network (LAN) or to a network management station over a 
telecommunication network. 

A final remark is that embodiments of the present invention are 
described above in terms of functional blocks. From the functional 
description of these blocks, given above, it will be apparent for a person 
15 skilled in the art of designing electronic devices how embodiments of these 
blocks can be manufactured with well-known electronic components. A 
detailed architecture of the contents of the functional blocks hence is not 
given. 

While the principles of the invention have been described above in 
20 connection with specific apparatus, it is to be clearly understood that this 
description is made only by way of example and not as a limitation on the 
scope of the invention, as defined in the appended claims. 




